Ground State Energy Fluctuations of a System Coupled to a Bath 
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It is often argued that a small non-degenerate quantum system coupled to a bath has a fixed 
energy in its ground state since a fluctuation in energy would require an energy supply from the 
bath. We consider a simple model of a harmonic oscillator (the system) coupled to a linear string 
and determine the mean squared energy fluctuations. We also analyze the two time correlator of 
the energy and discuss its behavior for a finite string. 



The coupling of two subsystems is a central problem 
in many areas of physics. The interaction of atoms with 
the radiation field haS|-been central in the development 
of quantum mechanicsEJ. The division into subsystems 
often depends on the questions asked: For instance de- 
phasing in mesoscopic systems investigates the quantum 
motion of a single electron in a conductor viewing all 
other conduction electrons as the bathH. fn this work 
we are interested in the energy fluctuations in a test sys- 
tem coupled to another system. Of particular interest 
is the zero-temperature limit. It is often argued that in 
the zero-tempieEature limit two subsystems can not ex- 
change energyOa. The argument is based on the assump- 
tion that both subsystems are in their separate ground 
state (which they would assume in the absence of any 
coupling). Since both the test system and the bath are 
in their ground sta|tjft, neither of the two can supply an 
energy to the otherlju. Below we consider a simple sys- 
tem of an oscillator (the test system) coupled to a linear 
string (the bath) and investigate the energy fluctuations 
of the test system. This is an exactly solvable model 
which demonstrates the existence of energy fluctuations 
in the zero-temperature limit. These fluctuations are a 
consequence of the flnite coupling energy between the 
test system and the bath. 

Often the interaction of two subsystems is treated in 
terms of the Einstein coefficients for the absorption and 
the spontaneous and stimulated emission of the test sys- 
tem. In such an approach rates of transitions of the test 
system to change energy form one eigenstate En to an- 
other eigenstate Em are investigated. These are inelastic 
transition rates which vanish when both systems are in 
the ground state. The resulting probabilities are used 
in a classical stochastic equation for the probabilities of 
occupation of the test system. Only probabilities en- 
ter in the description of the reservoir system interaction 
and this approach again yields no energy fluctuations in 
the ground state^ Such a description is challenged both 
in Laser physicso and in the morg recent discussions of 
macroscopic quantum cohcrenced'B. 

The role of zero-point fluctuations in mesoscopic con- 
ductors is a hotly debated issue and we cite here only 
Ref. ^ as an entry to the literature. The point of view 
taken here has been applied to investigate the persistent 



current of a small mesoscopic loop with a|-Quantum dot 
capacitively coupled to a transmission linea. The persis- 
tent current is a measure of the quantum coherence of 
the ground state. It was found that the ground state un- 
dergoes a crossover form a state with a well defined per- 
sistent current much larger than its mean square fiuctua- 
tions to a state in which the magnitude of the persistent 
current is much smaller than its mean-square fiuctuations 
as the coupling to the bath increaseaS. 

Vacuum fiuctuations can have important effects on the 
system considered. The Lamb shift is a widely appreci- 
ated example. Another example is the Casimir effect. 
In both of these examples the effect of the vacuum can 
be thought of in terms of a rcnormalizatifln of the orig- 
inal parameters characterizing the systenJj. In contrast, 
the energy fiuctuations which we discuss here can not be 
obtained from the uncoupled system simply by renormal- 
izing its parameters: After all, a harmonic oscillator in 
its ground state does not fiuctuate. 

The system we consider is a harmonic oscillator with 
mass mo and frequency uq described by the energy op- 
erator 



Ho 



2mo 



P ^ TT^I ■ 



(1) 



The oscillator is coupled to a harmonic bath, which leads 
to a dissipation linear in its velocityQ with a friction con- 
stant rj. The expectation value of the mean energy of the 
test oscillator in the ground state of the system is to the 
linear order in the friction constant ry 
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and the mean squared energy (SHq) = {{Ho — {Hq))'^) is 
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where uJc is a (high-frequency) Debye cut-off of the spec- 
trum of the bath. Thus both the mean energy and the 
mean squared fiuctuations in energy increase (for weak 
coupling) linearly with the coupling constant 77. Note 
that the mean squared fiuctuations are proportional to 
. The result for the mean energy is well knownB: 
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In contrast, the result for the energy fluctuations seems 
novel. Below we extend Eq. to all orders in rj and 
investigate also the two time correlator of the energy (a 
fourth order correlation). Eqs. (H) and (||) are valid for 
the infinitely long string. To emphasize that the overall 
energy (system plus coupling energy plus bath) is con- 
served despite the energy fluctuations in the test system 
alone, we will later also consider a string with a finite 
number of particles only. 

The mean squared energy fluctuations can be deter- 
mined from the coordinate correlation function of the 
oscillator J3 which is conveniently obtained using the clas- 
sical response of the oscillator and the fluctuation - dis- 
sipation relation at a flnite temperature T 
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for luq < ri/2. The energy fluctuations are determined by 
the fourth order correlations of the momentum and coor- 
dinate. The momentum is related to the coordinate via 
p = moq and thus the momentum correlations can be re- 
duced to time-derivatives of coordinate correlations. We 
next make use of the fact that the fluctuations are Gaus- 
sian and thus fourth order correlations can be expressed 
as sums of products of second order correlations. Thus 
we obtain the energy correlator in the form 
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and then passing to the limit T = 0. Introducing $7+ = 
(wg - ?7V4)'/', for uJn > r?/2 and = {rf/4 - ^2)1/2 
for loq < ri/2- Eq. (||) gives at < = the mean squared 
displacements 
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For the mean squared fluctuations we need the time 
derivatives of the correlator at t = 0. Evaluating the 
resulting integrals yields, 



UJQ 



+ -mo77(cJo - ry2/2) In 



LUq 



„2\2 



-'0' 



(8) 



where (q^) is given by Eq. (0). Expanding Eq. 
terms linear in rj gives Eq. (Q). The full result Eq. 
shown in Fig. |^ together with the mean energy 
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of the test oscillator. 

So far we considered the system consisting of the oscil- 
lator and the bath in a thermodynamic equilibrium at a 
given temperature T, which was then set equal to zero. 
It is also instructive to trace the origin of the energy 
fluctuations by using a microscopic model of the bath 
and purely quantum-mechanical considerations without 
involving any thermodynamic relations. 

It is well known that a harmonic bath caia, bejiiojieled 
by a semi-infinite string of identical particleaa'Bll3il3 with 
mass rrih and frequency LUh- The electrical analog of this 
system consists of an LC-oscillator coupled to a trans- 
mission lineEfl. In the continuous limit where m/i —> 
and ujh ^ oo, the string leads to dissipation that is lin- 
ear in the velocity of the oscillator with friction constant 
r] = (mh/TriQjujh- In this limit, the interaction with the 
string does not shift the resonance frequency of the test 
oscillator and results only in a finite damping. The role 
of the high-frequency cutoff ojc is played now by ojh- 



To address a system with a well defined ground-state 
energy, we consider a test system attached to a chain of 
a finite length consisting of N oscillators. The Hamilto- 
nian of a system coupled to a linear harmonic string is 
conveniently split into the following three parts 
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FIG. 1. The mean energy (Ho) (broken line) and the 
mean-squared energy fluctuations {SHq) in the ground state 
as a function of the couphng constant rj in units with h = 1. 
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A transition from a finite to the infinite system is most 
conveniently traced by watching the behavior of a two- 
time correlator of energy. We introduce the exact classi- 
cal normal modes of the total Hamiltonian Ti. and calcu- 
late the correlator (|^) . Note that this is not the quantum- 
mechanical diagonalization of Ti. in the basis of its eigen- 
states because the number of these states is infinitely 
large even for one oscillator. Instead, we perform a linear 
transformation of coordinates Xa with an7V-|-lxA^-|-l 
matrix U that will simultaneously bring the quadratic 
forms 
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describing the kinetic and potential energy of the system 
to a diagonal form. To this end, we introduce the classi- 
cal normal modes of the system ipk, which obey the set 
of equations 
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Here ^ — {mh/moy^^ is the ratio of the mass of the bath 
oscillators and the mass of the test oscillator. We search 
the eigenvectors in the form 



ipk{a >l)=Ak cos(Afca) + Bk sin(Afca). 
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On applying the transformation ( |l7| ) to the Hamiltonian 
TYl, it assumes the form 
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The quantum-mechanical coordinate of the test oscillator 
q can be written in terms of the normal coordinates ^fc 
of the transformed Hamiltonian Id dl8) 
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and the latter can be presented in the form 
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where cik and are time-independent annihilation and 
creation operators with the standard commutation rules. 

Now one easily obtains from Eq. (|^) an expression for 
the correlation function of energy fluctuations 
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FIG. 2. Symmetrized correlator of the energy fluctuations versus time for uio/coh = 1, rrih/mo — 0.1, and different numbers 
of particles A'^ in a finite string in units with h — 1. The left arrows indicate the value of the correlator at f = 0. 



The amplitudes Ak and eigenfrequencies ujk were nu- 
merically calculated and equation was evaluated for 
LUo/ujh = 1, to/j/too = 0.1, and different values of N. The 
results are shown in Fig. |[ 

Since the string is finite, we now have a recurrence phe- 
nomenon. As shown in Fig. |2[ the correlation function 
initially oscillates and its envelope decreases rapidly (as 
it would for the infinite chain), but it rebuilts after a time 
it takes a signal to travel along the string and to return 
to the oscillator after reflection at the opposite end. The 
built-up time is different from the decay time of the corre- 
lation and in detail the reconstituted correlation is clearly 
not the same as the initial correlation. This behavior is 
a consequence of the fact that the eigenfrequencies Uk of 
our problem are only nearly exact multiples of the low- 
est eigenfrequency. Even in the limit of a continuous yet 
finite string {nih —> 0, LUh ^ oo), the equidistant levels 
of the system are disturbed near the test-oscillator fre- 
quency uiQ. The fact that the eigenfrequencies are not 
exact multiples leads to the more complex phenomenon 
depicted in Fig. ^. The above results suggest that a tran- 
sition to the thermodynamic limit is possible if the travel 
time of a perturbation through the system is longer than 
the duration of an experiment. 

The exact diagonalization demonstrates that the en- 



ergy of a small system coupled to a bath fluctuates as a 
function of time despite the fact that the overall energy 
of the system is fixed. This is a consequence of the finite 
coupling energy between the system and the bath. As a 
result the test oscillator is not in a normal mode of the 
total system but as shown by Eq. (|l9| ) in a superposition 
of the true normal modes of the total system. It is for 
this reason that the simple argumentLiQ mentioned in the 
introduction fails. 

The results presented in this work demonstrate that 
even at zero temperatures vacuum fluctuations of the 
bath give rise to a ground state that exhibits a-jion- 
trivial dynamics. Textbook statistical mechanicsEj as- 
sumes that the coupling energy between the test system 
and the bath can be neglected. In contrast, here the cou- 
pling energy is essential. The effect discussed here can 
not simply be absorbed into renormalized parameters of 
the test system. We believe that these observations are 
crucial in understanding the zero-temperature properties 
of systems and that they are applicable not only to the 
particular system investigated here but are of a very gen- 
eral nature. 

This work was supported by the Swiss National Sci- 
ence Foundation and the Russian Foundation for Basic 
Research, grant 01-02-17220. 



4 



^ C. W. Gardiner and P. Zoller, "Quantum Noise" , (Springer, 

Heidelberg, 2000). 
^ B. L. Altshuler, A. G. Aronov and, D. Khmelnitskii, J. 

Phys. C 15, 7367 (1982) 
^ I. L. Aleiner, B. L. Altshuler, and M. E. Gershenson, Waves 

in Random Media 9, 201 (1999). 

U. Gavish, Y. Lcvinson, and Y. Imry, Phys. Rev. B 62, 

R10637 (2000). 
^ A. J. Leggett, et al.. Rev. Mod. Phys. 59, 1 (1987). 
® U. Weiss, Quantum Dissipative Systems, (Word Scientific, 

2000) 

^D. Natelson, et al., Phys. Rev. Lett. 86, 1821 (2001); J. 
J. Lin and L. Y. Kao, J. Phys. Condens. Matter 13, L119 



(2001). 

® P. Cedraschi, V. V. Ponomaronko, and M. Biittikcr, Phys. 
Rev. Lett. 84, 346 (2000); P. Cedraschi and M. Biittiker, 
Annals of Physics, 289, 1 - 23 (2001); Phys. Rev. B 63, 
165312 (2001). 

® P. W. Miloni, The Quantum Vacuum, (Academic Press, 

Boston, 1994). 
1° H. Lamb, Proc. London Math. Soc. 53, 208 (1900). 
" R. J. Rubin, J. Math. Phys. 1, 309 (1960); ibid 2, 373 
(1961). 

12 H. Grabert, U. Weiss and P. Talkner, Z. Phys. B55, 87 
(1984). 

" G. W. Ford, J. T. Lewis, and R. F. O'Connell, Phys. Rev. 

A 37, 4419 (1988). 
" B. Yurke and J. S. Denker, Phys. Rev. A 29, 1419 (1984). 

R. Becker, Theorie der Warme, (Springer Verlag, Berlin, 

1966). 



5 



